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3^P-NMR and ^SR Studies of Filled Skutterudite Compound SmFe4Pi2: 
Evidence for Heavy Fermion Behavior with Ferromagnetic Ground State 
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The ^^P-NMR (nuclear magnetic resonance) and /xSR (muon spin relaxation) measurements 
on the filled skutterudite system SmFe4Pi2 have been carried out. The temperature T depen- 
dence of the ^'^P-NMR spectra indicates the existence of the crystalline electric field effect 
splitting of the Sm'^"'" (J = 5/2) multiplet into a ground state and an excited state of about 
70 K. The spin- lattice relaxation rate 1/Ti shows the typical behavior of the Kondo system, 
i.e., l/Ti is nearly T independent above 30 K, and varies in proportion to T (the Korringa 
behavior, l/Ti oc T) between 7.5 K and 30 K. The T dependence deviated from the Korringa 
behavior below 7 K, which is independent of T in the applied magnetic field of 1 kOe, and 
suppressed strongly in higher fields. The behavior is explained as l/Ti is determined by ferro- 
magnetic fiuctuations of the uncovered Sm^^ magnetic moments by conduction electrons. The 
/iSR measurements in zero field show the appearance of a static internal field associated with 
the ferromagnetic order below 1.6 K. 

KEYWORDS: SmFe4Pi2, filled skutterudite, ferromagnetism, heavy fermion system, ^^P-NMR (nuclear 
magnetic resonance) , ^SR (muon spin relaxation) 
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1. Introduction 

Filled skutterudite compounds RT4X12 {R: rare earth, 
actinide; T: Fe, Ru, Os; X: P, As, Sb) crystallize in a 
body-centered cubic (BCC) structure of the space group 
Im 3 (No. 204).^ These compounds have recently at- 
tracted much attention as improved thermoelectric ma- 
terials and for their wide variety of electrical and mag- 
netic properties, such as the superconductivity, the heavy 
fermion (HF) system, the multipolar order, the magnetic 
order, the metal-insulator (M-I) transition and the non- 
Fermi-liquid behavior. 

Among them, SmRu4Pi2 was reported to exhibit an 
M-I transition at Tmi of 16.5 K. Recently, it was demon- 
strated that the M-I transition accompanies the magnetic 
octupolar order and that there is the additional antifer- 
romagnetic (AFM) order below 15 K, respectively.''"^^ 
The SmOs4Pi2 was reported to show the simple AFM 
order below 4.6 K.^^." The SmFe4Pi2 and SmOs4Sbi2 
were reported to exhibit the heavy fermion behavior with 
ferromagnetic (FM) ground states. ^'^"^^ 

In SmFe4Pi2, susceptibility and specific heat measure- 
ments show that the Sm ions are in a trivalent state 
and that the system has an FM transition at the Curie 
temperature Tq of 1.6 K.^^'^^ The temperature T de- 
pendence of the electric resistivity exhibits the Kondo- 
lattice behavior (Kondo temperature Tk = 30 K), and 
the electronic specific heat coefficient attains the value of 
370 mJ/mole K^.^^' in a BCC structure with the point 
group of T/i symmetry, the Sm'^^" (J — 5/2) multiplet is 
split into a doublet T<^ and a quartet Pgr by the crys- 
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talline electric field (CEF). The specific heat data also 
suggests the CEF effect splitting of the ground state and 
the excited state separated by around 70 K.^^ However, 
the Kondo effect brings the uncertainty in the determi- 
nation of the ground and the excited states, i.e., we do 
not know whether the ground state is Ps or Pg7. Indeed, 
the magnetic entropy associated with the FM order is 
0.17i?ln 2, which is much smaller than the value expected 
for the case that the ground state is P5. 

In this paper, we report the microscopic magnetic 
properties of SmFe4Pi2 obtained by '^^P-NMR (nuclear 
magnetic resonance) and /iSR (muon spin relaxation) 
measurements. The results of the "^^P-NMR Knight shift 
K and the spin-lattice relaxation rate l/Ti measure- 
ments are discussed on the basis of the Kondo screened 
magnetic moment model. The /iSR is ideal for the in- 
vestigation of the magnetic properties in zero external 
field. 

2. Experimental 

The single-phase polycrystalline SmFc4Pi2 has been 
synthesized by using the high temperature and high pres- 
sure method. ^'^ The sample was crushed into powder for 
the experiments. 

The •^^P-NMR has been performed by using phase- 
coherent pulsed NMR spectrometers and superconduct- 
ing magnets. The '^"'^P NMR spectra were measured by 
the convolution of Fourier transform signals of the spin 
echo which were obtained at 20 kHz interval. The 1/Ti 
was measured with the saturation recovery method in ap- 
plied magnetic fields of 0.921-33.4 kOe. The nuclear mag- 
netization recovery curve was fitted by a simple exponen- 
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tial function as expected for the nuclear spin / = 1/2 of 
^^P nucleus, which allows us to determine a unique 1/Ti 
value at each temperature and field. 

The /LtSR has been performed by implanting pulses of 
positive surface muons at the RIKEN-RAL Muon Facil- 
ity in the UK. The direction of the initial muon spin is 
parallel to the beam line where forward and backward 
counters were located on the upstream and the down- 
stream sides. The asymmetry parameter A{t) was de- 
fined as A{t) = [F{t) - aB{t)]/\F{t) + aB{t)] - Abg, 
where F{t) and B{t) are the total muon events counted 
by the forward and the backward counters at time t, re- 
spectively. The Abg is the background, and a is the cali- 
bration factor reflecting the relative counting efficiencies 
of both counters. In zero field (ZF-) experiments, stray 
fields at the sample position were compensated within 
0.03 Oe by using correction coils, which is small enough 
for our ZF-/xSR measurement. 

3. Results and Discussion 

Figure 1 shows the ^^P-NMR spectra at several tem- 
peratures in the applied field of 94.063 kOe. The line 
shape of the spectra at high temperatures shows a pow- 
der pattern with a nearly uniaxial Knight shift distri- 
bution for nuclear spin I — 1/2. However, the line shape 
gradually changes and deviates from the uniaxial powder 
pattern with decreasing temperature. The change of the 
line shape becomes remarkable below 70 K. The spectra 
were analyzed by taking into account of the powder aver- 
age of K anisotropy and the excess Gaussian broadening. 
The K is expressed as 

K = Kiso + ^filScos^ 0-1)- JCalsin^ 6'cos2(?i), (1) 

where the Kiso represents an isotropic term, Ki repre- 
sents a uniaxial term, and K2 represents the deviation 
from the uniaxial symmetry. 

The T dependence of Xiso, —Ki and K2 are shown in 
Figs. 2 and 3, respectively. The T dependence of Kiso is 
remarkably different from that of the bulk susceptibility 
X- The Kiso rapidly increases with increasing tempera- 
ture, and has a broad maximum around 120 K and then 
decreases slightly at higher temperatures. To clarify the 
T dependence, Kiso is plotted against x, as seen in Fig. 4. 
The K-x plot shows the T dependence of the transferred 
hyperfine coupling as 

where iVA and /iB are the Avogadro number and the 
Bohr magneton, respectively. It is noticed that the K- 
X plot is linear in the T ranges of 25-120 K and above 
about 150 K. The hyperfine coupling constants Ahf were 
estimated about —0.71 kOe/nB below 120 K and about 
-1-0.74 kOe//XB above 150 K from the K-x plot by adopt- 
ing the formula (2). Considering Tk = 30 K in this sys- 
tem, the change in the hyperfine coupling around 120 K 
is reasonably attributed to the CEF splitting. The change 
of a hyperfine field associated with a CEF effect was re- 
ported in other system. 

The Ki is nearly T independent in the whole T region 
with the value of about —0.015 %. The K2 decreases 
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Fig. 1. ^^P-NMR spectra at several temperatures in the applied 
magnetic field of 94.063 kOe. The broken line shows the reso- 
nance frequency of 162.14 MHz at K = Q estimated from the 
nuclear gyromagnetic ratio of the ^^P nucleus. 



with increasing temperature below about 70 K, and de- 
creases very gradually at higher temperatures (0.0059 % 
at 280 K). One origin of the anisotropic hyperfine cou- 
pling is the the direct dipolar coupling of Sm'^+ mag- 
netic moments and ^^P nuclear spins. The lattice sum 
calculation predicts nearly uniaxial Knight shift to be 
Ki ~ 0.014 % and K2 ~ 0.0017 % at 280 K. These 
results are shown in Table I including the experimen- 
tal results. Though the uniaxial anisotropy is observed 
in our system at high temperatures, Ki has the oppo- 
site sign. This result indicates that the hyperfine mech- 
anism can not be explained by only the direct dipolar 
coupling even at sufficiently high temperatures (280 K). 
The nearly same magnitude of K2 at high temperatures 
between the experimental and the calculated values is ac- 
cidental. This suggests that another hyperfine field mech- 
anism contributes to the anisotropic Knight shift. The 
hybridization of the Sm-4f orbital with the P-3p orbital 
makes spin polarizations at P-3p spins, which induces 
the anisotropic field via the dipolar coupling of the P- 
3p spins and ^^P nuclei. The sum of these contributions 
would explain the observed values. The change of the 
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Fig. 2. Temperature dependences of Ki^ 
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Fig. 3. Temperature dependence of —Ki (closed circles) and K2 
(open circles). 



spectra below 70 K is represented by the increase of K2 
which reflects the CEF splitting of about 70 K.^* The 
change of a Sm-4f orbital shape occurs associated with 
the thermal excitation of the CEF levels and modifies the 
hybridization of the Sm-4f and the ligand P orbitals. For 
the evaluatation of the ligand spin polarization, a com- 
plete solution of the hybridization process is required. 
The analysis of the spectra is an open problem for future 
works. 

Figure 5 shows the T dependence of 1/Ti in several ap- 
plied fields. The Ti tends to become long with decreasing 
temperature. Above 30 K, the weak T dependence was 
observed. This is attributed to the relaxation by the 4f 
localized moments of the Sm ions. We evaluated the con- 
tribution of the conduction electrons by 1 /Ti of the non- 
magnetic LaFe4Pi2, which is much smaller than the value 



Fig. 4. Bulk susceptibility dependence of Kiso {K-x plot). The 
solid lines show the linear function of x expressed by the formula 
(2). Some temperatures are also shown corresponding to x- 



Table I. Experimental and the dipolar calculated values of Ki 
and at the temperature of 280 K. 





Ki (%) 


K2 (%) 


Experimental 


-0.015 


0.0059 


Calculated 


0.014 


0.0017 



observed in SmFe4Pi2.^" After subtracting this contribu- 
tion, the corrected 1/Ti has much weaker T dependence, 
which is shown as "(l/ri)4f" in Fig. 5. We observed the 
Fermi liquid (the Korringa behavior, TiT ~ constant) in 
the temperature range of 7.5-30 K, which indicates the 
major part of the Sm 4f-electrons to be itinerant by the 
hybridization of the conduction electrons below Tk- 

The 1/Ti at 0.921 kOe is nearly T independent above 
Tc, and rapidly decreases at lower temperatures. The 
rapid decrease suggests the existence of a phase tran- 
sition in ZF. The critical slowing-down phenomenon of 
1/Ti might be suppressed even by the small external 
field. We estimated 1/Ti on the presumption that the 
Kondo screened Sm^+ moments ^0.2/13 coupled by the 
ferromagnetic exchange interactions.^^ The 1/Ti is ex- 
pressed as 
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2zj{j + iy 

J is the quantity of the angular momentum of the Sm 
ions, Z is the number of the nearest neighbor Sm ions, h 
is the Planck constant divided by 27r, 7 is the nuclear gy- 
romagnetic ratio of the ^^P nucleus, ffint is the internal 
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Fig. 6. ZF-/iSR time spectra at several temperatures. The solid 
lines show the best-fit results by the formula (6). 



Fig. 5. Temperature dependence of 1/Ti in several applied mag- 
netic fields. The 1/Ti of nonmagnetic LaFe4Pi2 is also shown 
as a reference.'^" "(1/Ti)4f" is obtained by subtracting 1/Ti of 
LaFe4Pi2 from that of SmFe4Pi2. 



field made by the 4f localized moments of the Sm ions at 
the ^^P spins, /U is the saturated moment at 1.4 K, J^xc 
is the ferromagnetic exchange interaction energy, and fee 
is the Bolzmann constant. The estimated value of 1/Ti 
is about 11.6 1/s, which is close to the experimental val- 
ues. This result stiggcsts that the origin of the ferromag- 
netic ground state would be due to the uncovered 4f local 
moments of the Sm^+ by the conduction electrons. The 
value of 1 /Ti at low T decreases rapidly with increasing 
field. The coupling of the uniform ferromagnetic fluctu- 
ations and external field strongly suppresses 1/Ti. 

It is important to carry out ZF-measurements to iden- 
tify the anomaly around Tq as a phase transition. How- 
ever, it is impossible to perform ^-"^P-NMR in ZF, since 
the ^"'^P nucleus (7=1/2) has no quadrupole moment. 
Hence, we carried out ZF-/xSR to confirm the phase tran- 
sition. 

Figure 6 shows the ZF-//,SR time spectra at several 
temperatures. Below Tc, A{t) rapidly decreases, and the 
muon-spin precession was observed as seen clearly in 
Fig 7. This indicates that the system has a magnetically 
ordered ground state. The internal field of about 650 Oe 
at the muon site was estimated from the precession fre- 
quency by the formula (6). Above Tc, A{t) gradually 
decreases with t where muon-spins depolarize slowly by 
transferred magnetic fields from dynamically fluctuating 
Sm moments and nuclear moments. 

The spectra were analyzed by a multi-components 
function expressed by the formula 

A{t) = Aie-^i* + Ase"^^* -|- ^36-^=** cos(a;i -|- 9). (6) 




Fig. 7. ZF-/1SR time spectrum at the temperature of 0.31 K, 
which is obtained at a high time resolution. The solid line shows 
the best-fit results by the formula (6). 



The first and the second terms simply represent the com- 
ponents of the rapid and the slow depolarizations, re- 
spectively. The third term describes the precession com- 
ponent. The parameters of Ai, A2, A3 are the initial 
asymmetries at t = 0, the Ai, A2, A3 are the muon-spin 
depolarization rates, and oj, 9 are the frequency and the 
phase of the precession, respectively. The spectra above 
1 K which have no precession were analyzed without the 
third term. Figures 8 and 9 show the T dependence of 
Ai, Ai and A2, A2 above 1 K obtained from the analysis, 
respectively. It is noted that the sum of Ax and A^ is the 
constant value of A{Q). The spectra above 3 K are almost 
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Fig. 8. Temperature dependence of Ai (closed circles, left axis) 
and Ai (open circles, right axis) above 1 K obtained from the 
best-fit results shown as the solid lines in Fig. 6. 
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Fig. 9. Temperature dependence of A2 (closed circles, left axis) 
and A2 (open circles, right axis) above 1 K obtained from the 
best-fit results shown as the solid lines in Fig. 6. 



dominated by only the second term as seen in Ai and 
A2 where the spectra almost follow a single exponential 
function. Hence, the absolute values of Ai (correspond- 
ing to A2) and Ai above 3 K are not so intrinsic, and the 
value of Ai should be idealy. Below 3 K, the first term 
(rapid depolalization) represents the relaxation by both 
static and dynamical components which develop below 
the temperature. Then, we forcus our attention on A2 of 
the second term to discuss the dynamical component. 

As seen in Fig. 9, clear anomalies were observed in 
both A2 and A2 around Tc- The rapid decrease of A2 
corresponds to the rapid increase of Ai. The appearance 
of the fast relaxation behavior suggests the development 



of a static internal field attributed to the existence of a 
magnetic order. Considering that the slow relaxation of 
the spectra is due to a dynamical component, the peak 
in A2 arises from the critical slowing down behavior of 
the Sm moments around Tq. The observation of the crit- 
ical slowing down behavior seems to contradict that of 
1/Ti by "^^P-NMR. One reason would be due to the time- 
window difference of each probe. The gyromagnetic ratio 
of the muon is two orders of the magnitude larger than 
that of ^^P nucleus. Hence, /iSR is much sensitive to 
higher frequency components of magnetic fluctuations, 
which makes the broad temperature dependence of A2. 
It is noteworthy that the behavior of A2 is not due to the 
inhomogeneity of the sample because the sample used in 
this measurement had been well characterized by other 
macroscopic measurements. 

Though the bahavior of A2 and A2 is obtained quali- 
tatively, this is consistent with the results of ^^P-NMR 
and other measurements.^^' Indeed, these behavior 
was also observed around T^i in SmRu4Pi2 where the 
magnetic order below Tmi has been confirmed by previ- 
ous /xSR measurements.^ The most important point of 
this ZF-^SR result of SmFe4Pi2 is the observation of the 
muon-spin precession below Tc. These results conclude 
that the magnetic order occurs below Tc, and support 
the ^^P-NMR results which could not be obtained in ZF. 

The FM fluctuations around Tc in finite fields and the 
existence of a magnetic phase transition at Tc in ZF were 
observed by means of ^^P-NMR and fiSR, respectively. 
Therefore, an FM phase transition at Tc in SmFe4Pi2 
has been confirmed. 

4. Conclusions 

The 3ip-NMR in external fields and the ZF-/iSR stud- 
ies on SmFe4Pi2 have been carried out. These results 
show the existence of the CEF effect, the HF behavior 
below Tk of 30 K and the FM order below Tc of 1.6 K. 
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